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Abstract 

Using supersymmetric grand unified theories, we have recently in- 
vented a framework which allows the prediction of three quark masses, 
two of the parameters of the Kobayashi-Maskawa matrix and tan /3, the 
ratio of the two electroweak vacuum expectation values. These predic- 
tions are used to calculate e and e' in the kaon system, the mass mixing 
in the i?^ and systems, and the size of CP asymmetries in the decays 
of neutral B mesons to explicit final states of given CP. 



*This work was supported in part by the Director, OfRce of Energy Research., Office of 
High Energy and Nuclear Physics, Division of High Energy Physics of the U.S. Department 
of Energy under Contracts DE-AC03-76SF00098 and DOE-ER-01545-573 and in part by the 
National Science Foundation under grants PHY90-21139 and PHY86-12280. 



Disclaimer 



This document was prepared as an account of work sponsored by the United States Gov- 
ernment. Neither the United States Government nor any agency tiiereof, nor The Regents 
of the University of California, nor any of their employees, makes any warranty, express or 
implied, or assumes any legal liability or responsibility for the accuracy, completeness, or use- 
fulness of any information, apparatus, product, or process disclosed, or represents that its use 
would not infringe privately owned rights. Reference herein to any specific commercial prod- 
ucts process, or service by its trade name, trademark, manufacturer, or otherwise, does not 
necessarily constitute or imply its endorsement, recommendation, or favoring by the United 
States Government or any agency thereof, or The Regents of the University of California. The 
views and opinions of authors expressed herein do not necessarily state or reflect those of the 
United States Government or any agency thereof of The Regents of the University of California 
and shall not be used for advertising or product endorsement purposes. 



Lawrence Berkeley Laboratory is an equal opportunity employer. 



ii 



In previous papers [1] we have invented a predictive framework for quark 
and lepton masses and mixings based on the Georgi-Jarlskog ansatz [2] for the 
form of the mass matrices in supersymmetric grand unified theories [3]. In this 
paper we use this scheme to make predictions for parameters of the neutral K 
and B meson systems. In particular we show that the CP asymmetries in neutral 
B meson decays are large and will provide a powerful test of the scheme. We 
begin by reviewing the predictions for quark masses and mixings. 

The top mass is predicted to be heavy 

^..f rrib \f \ /.053\ Vl-46\ /1.84\ 

where rii is the QCD enhancement of a quark mass scaled from vtit to vtii. Per- 
turbativity of the top Yukawa coupling also requires that vtit <187 GeV. In this 
paper the central values of r]i quoted correspond to a complete two loop QCD 
calculation with UgiMz) = -109, whereas in reference 1 an approximate two loop 
result was given. 

The value of = .109 comes from a 1 loop analysis of the unification of 
gauge couplings which for simplicity ignored threshold corrections at the grand 
unified and supersymmetry breaking scales [1]. However these threshold correc- 
tions will be present at some level in all grand unified models [4], and would not 
have to be very large for our predicted value of the QCD coupling to range over 
all values allowed by the LEP data: .115± .008. Whatever the threshold correc- 
tions are, they must give an acceptable value for the strong coupling. Hence it is 
important to consider the range of top masses allowed by the LEP range of a^. 
Larger values of lead to larger rji reducing rut. Increasing as{Mz) from .109 
to .123 decreases mt from 179 GeV to 150 GeV. Alternatively, as{Mz) = .123 
allows the top mass to be near the fixed point value of 187 GeV with Vcb = .047. 
The above numbers refer to the running mass parameter. The pole mass, which 
is to be compared with experiment, is 4.5% larger. 

The particular form of the quark mass matrices leads to an unusual form 
for the Kobayashi-Maskawa matrix 

/ Ci - SiS2e"*'^ Si + CiS2e~"*' S2S3\ 

V = -C1S2 - sie-'^ cie-'^ - S1S2 S3 (2) 
V siss -C1S3 e''*' ) 

where si = sin^i, etc, and we have set C2 = C3 = 1. f\ We do not lose any 
generality by choosing the phases of quark fields such that 61,62 and 6^ all lie 
in the first quadrant. We have predicted [1] 

si = .196 



*The angle 63 used in this paper corresponds to 63 — 64 used in reference 1. 
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S2 = -OSSx 



(3) 



where 

_ j iriu/md l.22GeV r]c " , . 

^ ~ V 0.6 rric 1.84 ^ 
while the input Vcb determines S3 which must be chosen quite large in view of (1). 
The ratio ms/md, which we predict to be 25.15, strongly prefers niu/md < .8, 
while the present value of iT^b/Kbl — •^2 prefers mu/md larger than .4. In all 
of our predictions the largest uncertainty lies in mu/md, which we will display 
through the parameter x- 

The angle (f) is determined by the requirement that \Vus\ — sin 9c 

04,=- (0.51 (1 ± 0.11) - o.i3x^) = o.sstil 

~ 0.92 1^1.15 - ^ (1 ± 0.22) j = 0.92;:^^ (5) 

In the first expressions the x dependence is shown explicitly, together with the 
uncertainty from the measured value of the Cabibbo angle sin^^c = 0.221 ± .003. 
Note that the 0(1%) uncertainties in 9c become greatly magnified in 0. For 
this reason we keep track of the 9c dependence in our predictions. For cos0 
the expression is exact, while for sin0 it is good to better than l%.The final 
numerical expressions correspond to the limits x^ = 1=f| and sin^c = .221±.003, 
which are used for all numerical predictions in this paper. Notice that is 
determined to be positive and the experimental data on Re e in the kaon system 
forces positive. Hence there is no quadrant ambiguity: choosing 9i^2,3 all 
in the first quadrant means that (j) is also in the first quadrant. The rephase 
invariant measure of CP violation [5] is given in our model by 

J = ImVudVtbV*t,V^*a = CiC2C3SiS2sls^ 

where 

fix) = X (i-i5 - ^ (1 ± 0-22)) = i;:i. (7) 

The scheme which leads to these predictions involves mass matrices at the uni- 
fication scale with seven unknown real parameters. Six of these are needed to 
describe the eigenvalues: m„ ^ ^ and md ^ ^ m^, while the 
seventh is the CP violating phase. Hence a more predictive theory, having fewer 
than seven input parameters, must either relate the up mass matrix to that of 
the down, or must have an intrinsic understanding of the family mass hierarchy. 
Without solving these problems the most predictive possible theory will involve 
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seven Yukawa parameters. Such a predictive scheme can only be obtained by 
relating the parameters of the lepton mass matrix to those of the down quark 
mass matrix. To our knowledge the only way of doing this while maintaining 
predictivity is to use the ansatz invented by Georgi and Jarlskog [2] . The crucial 
point about our scheme for fermion masses is that it is the unique scheme which 
incorporates the GUT scale mass relations = ttIt, rUg = and rrid — SrUg 

with seven or less Yukawa parameters and completely independent up and down 
quark matrices. The factors of three result from there being three quark colors. 

It is because of this uniqueness that the detailed confrontation of this model 
with experiment is important. If the model is excluded, for example by improv- 
ing measurements of mt,Vcb or Vub/Vcb, then the whole approach of searching 
for a maximally predictive grand unified scheme may well be incorrect. Alter- 
natively it may mean that there is a very predictive scheme, but it involves 
relations between the up and down matrices in an important way. It is im- 
portant to calculate the observable parameters of the neutral K and B meson 
systems as accurately as possible, so that future experiments and lattice gauge 
theory calculations will allow precision tests of this scheme. 

Our form for the Kobayashi-Maskawa matrix, equation (2), is very unfamil- 
iar and so we give predictions for the parameters that appear in the Wolfenstein 
form of the matrix [6]. The Wolfenstein form is an approximate form for the 
matrix which is unitary only to order A"^, where A = sin^c = {Vusl- To order A^ 
we find the matrix can be written as: 



/ l-f-^ A AX\a-tf3)\ 

-A 1-^-^-^ AX^ 

\A\^-A\\a + ip) -A\^ + ^ J 



where A is defined by Vcb and a, Phy Vub- The reason that we prefer to work with 
the matrix at O(A^) is that for us the Vub entry numerically really is O(A^). Thus 
we have A, a, (3 = 0(1). Notice that the order A^ contributions to Vud-, Ks, 
and Vtb can be dropped unless an accuracy of greater than 2 1/2% is required. 
This means that F^^^ is actually the same as V^^\ the usual Wolfenstein form 
with p — \a and rj — X(3. Hence we use the usual Wolfenstein parameters 
A, p, rj. We find 

A = 1.09^ 
.053 

p = ^ = 0.12(1 - 0.25x^) 

,7 =^ = 0.22/(x). (8) 

The leading dependence on the uncertain quantity mu/m^ is shown explic- 
itly, through the parameter x- The uncertainties in p and rj coming from sin6'c 
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are less than 10% and 4% respectively. These are not shown as the Wolfen- 
stein approximation is itself only good to about 20%. Notice that we can write 
p + ir] — re**^, where 

r = ^ = 0.24X. (9) 
Requiring unitarity for the imaginary part of V to O(A^) [6] one deduces 

J A'X% = 3.1 X 10-5 (^)' /(X). (10) 

This agrees well with the exact result of equations 5, since (10) is expected to 
have 0(A) ~ 0(20%) corrections. 

Beneath the scale of grand unification our effective theory is just that of 
the minimal supersymmetric standard model (MSSM). Hence in the rest of this 
paper we wish to give the predictions for e, e', x^, Xg and the CP violating an- 
gles a, 7 in the minimal supersymmetric standard model, with the Kobayashi 
Maskawa matrix given by equation 2, and with our predicted values for quark 
masses. 

In reference 7 it will be shown that in the MSSM the supersymmetric con- 
tributions to e^X(i-Xs and to the CP violating angles and 7 in B meson 
decay are small. Here we will simply give a simplified discussion of why these 
contributions are negligable for quark masses and mixings of interest to us. In 
the standard model with a heavy top quark the quantities e, Xd and Xs arc dom- 
inated by box diagrams with two internal top quarks. The amplitude of these 
standard model box diagrams can be written as 

Bij^AsM{VuV:^'' (11) 

where i, j = d,s,b label the relevant external mass cigcnstatc quark flavors of 
the diagram and the dependence on the Kobayashi-Maskawa matrix elements 
is shown explicitly. The leading supersymmetric contributions to these three 
quantities come from box diagrams with internal squarks and gluinos. In this 
case the flavor changes occur through off-diagonal squark masses: Mfj. The 
amplitudes for these box diagrams can be written as 

Blj ^ Amssm (Mfj)' (12) 

where again only the relevant flavor structure has been shown explicitly. It has 
been assumed that squarks of flavor i and j are degenerate. 

In the MSSM the squarks are all taken to be degenerate at the grand unifled 
scale. The squark mass matrices evolve according to renormalization group 
equations which generate non-degeneracies and flavor-changing entries. When 
all quarks are light a very simple approximation for the flavor-changing entries 
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of the mass matrices results [8] . Since we predict a top Yukawa coupling close to 
unity, this is not good enough for our purposes. We use the analytic solutions of 
the renormalization group equations valid to one loop order in the top Yukawa 
couphng, but with other Yukawa couplings neglected [9]. This is the same 
approximation used to obtain our quark mass and mixing predictions [1] and is 
sufficient providing tan f3 is not so large as to make the bottom Yukawa coupling 
large. In this approximation only SU(2) doublet squarks have flavor changing 
masses. We are able to find a very convenient approximation for the induced 
flavor changing mass squared matrix elements for the down type doublet squarks 



where M is the mass of the (nearly) degenerate squarks and ^ is the ratio of the 
gluino to squark mass at the grand unified scale. For values of the top Yukawa 
coupling consistent with the prediction of equation 1, this approximation is good 
to better than a factor of two. The exact result has only shght sensitivity to the 
trilinear scalar coupling A of the MSSM, which we have neglected. 

Comparing the standard model box amplitude (equation 11) to that of 
the MSSM superbox amplitude (equations 12 and 13) it is apparent that the 
dependence on the Kobayashi-Maskawa matrix elements is identical. Hence the 
ratio of box diagrams is independent of the external flavors i,j 



from the momentum integral of the superbox diagram [8] and varies from 1(1) 
= 1/30 to 1(0) = 1/3. If the squarks are taken light (eg 150 GeV) in an attempt 
to enhance the superbox amplitude, then x > 1 to avoid an unacceptably light 
gluino, resulting in / < 1/30. To increase I therefore requires an increase in 
M, but this rapidly decreases the importance of the superbox diagram. The ^ 
dependent factor in equation 14 is always less than unity. We conclude that the 
supersymmetric contributions to e, and x^ are unimportant in our scheme. 

We have shown that, in the MSSM with degenerate squarks, the superbox 
diagrams have the same Kobayashi-Maskawa phases as in the standard model 
box diagrams. This implies that the supersymmetric diagrams do not affect the 
CP asymmetry parameters a,f3 and 7 in B° decay, as is well known [10]. 

We have assumed that at the grand unified scale the squark mass matrices 
are proportional to the unit matrix. In supergravity theories this proportion- 
ality is expected only at the Planck scale. In general it is possible that large 
interactions of the quarks with superheavy fields could introduce large fiavor 




(13) 





where x = rrig/M and rrig is the gluino mass. The monotonic function I results 
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changing effects from renormalisation group scaling between Planck and grand 
scales [11]. We assume that this does not happen, as would be the case if the 
only large Yukawa coupling in the grand unified theory is that which generates 
the top quark mass. 

We now proceed to our predictions. Since the — Kg mass difference 
receives large long distance contributions we do not think it provides a useful 
test of our theory. On the other hand, all observed CP violation is described by 
the single parameter e, which is reliably calculated from short distance physics. 
To calculate e precisely we do not use the Wolfenstein form for the KM matrix. 
Instead we use a manifestly phase invariant formula for e in terms of J [12]. We 
find that the box diagram with internal top quarks dominates. Including a 20% 
contribution from the diagrams with one top and one charm we find 

The parameter Bk describes the large uncertainty in the matrix element of a 
four quark operator between kaon states. We use experiment for |e| and make a 
prediction for Bk'- 




(16) 

where we have used equations 1 and 6 for nit and J. The first uncertainty shown 
comes from the experimental value of |e| = (2.26 ± .02) x 10^^ while the second 
comes from sin^^c- This prediction for B^ is strikingly successful. Wc stress that 
mu/fad cannot vary too much in our theory: values larger than 0.6 arc strongly 
disfavored by the fact that m^/m^ is 25, while the present experimental values 
of Viib/Vcb disfavors niu/nid lower than 0.6. Allowing = 1 ± 1/3, sin^^c = 
0.221 ±.003 gives the range Bk = 0.31-0.57 (for rrih = 4.15 GeV and = 1.22 
GeV). This should be compared with recent lattice results B^ = 0.7 ± 0.2 in 
the quenched approximation [13]. 

In the standard model, predictions for e'/e are very uncertain because they 
depend sensitively on: i) various strong interaction matrix elements, ii) the value 
of Aqcd, iii) the value of the strange quark mass and iv) the value of the 
top quark mass. We use our central predictions for Aqqd and m^, and rely on 
the 1/N approximation for the QCD matrix elements [14]. Using the analytic 
expression given in reference 14 we are able to derive our prediction 

^ ~ 3.9 X 10"^ (2.7mt -^ + O.Sm^^ - 2.2mf) (^) (17) 

where is to be given in units of 179 GeV and is predicted in equation 1, and 
Bk is predicted in equation 16. This small result is not unexpected, given the 
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large value of rrit [15], and we stress that it is uncertain because we do not know 
how well to trust the 1/N matrix elements. 

The leading supersymmetric contribution to e'/e comes from a diagram 
with an internal gluino and an insertion of the flavor-changing squark mass of 
equation 13. We find that this superpenguin amplitude is small compared to 
the ordinary penguin: 

A(superpenguin) /'150GeV\^ / 1 + 3^^ \ ( a f ^ 

A{penguin) ~ ' V M ) 1^1 + 5.5^^ V ~ \l80GeV ) J 

(18) 

for the case of degenerate squarks and gluino of mass M. This is partly because 
the loop integral is numerically smaller, but is also because the ordinary penguin 
diagram is enhanced by an order of magnitude by a large Inm^ factor. Even 
though these supersymmetric contributions are negligable, the uncertainties in 
the QCD matrix elements still imply that e'/ e cannot be considered a precision 
test of our scheme. 

The dominant standard model contribution to mass mixing arises 

from the box diagram with internal top quarks. We find that for the B^: 

Using equation (1) for rrit, Vta — siSa and the experimental value for xa of .67±.10 
we predict 

■' \ rUb J \ J Vl.46/ Vl.84y .053 ^ ^ 

Wc note that if Xd,Tnt and Vtd are allowed to range over their experimentally 
allowed values, the prediction of the box diagram (equation 19) implies that 
VbJb will have to range over an order of magnitude. It is therefore a non-trivial 
success for our theory that it gives a prediction for \fBfB which is close to the 
quoted values for fB\fB- Our prediction should be compared with recent lattice 
results: fs = 205±40 MeV and v^/b = 220±40 MeV [16]. Our predictions for 
mt (1), Bk (16) and Jb (20) all depend on r]i which depend on ag- The numbers 
quoted are for the fairly low value of ag — .109. Threshold corrections at the 
grand unified scale could increase this, easily resulting in a 20% increase in r]br]c- 
This not only reduces the top mass, but gives improved agreement with lattice 
calculations for both Bk and /g. At any rate, once the top mass is measured 
our predictions for Bk and fs will be sharpened considerably. Our results for 
Bk and fs agree with those obtained reference 17. 

The standard model box diagram relates the mass mixing in B^ to that in 
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B'a by: 



r - w Pis p ; -^Ms^ ^2^^ 



3d/ X — ^JOd/ 

„2 /„2 



where we used our result for the ratio of Kobayashi-Maskawa factors: ci/sj — 
25 with neghgible uncertainty. If the ratio of B meson decay constants could be 
accurately calculated, and if large values of Xg (say 15 to 25) could be measured, 
then (21) provides a precision test of our theory. 

Finally we consider CP asymmetries which result when B'^ and B^ can 
decay to the same CP eigenstate / [18]. Unitarity of the KM matrix implies the 
1st and 3rd columns are orthogonal: VudV*/, + VcdYda + ^tdYth — 0- This can be 
represented as a triangle since the sum of three vectors is zero. Labelling the 
angles opposite these three vectors as a, and 7 respectively, one finds that 
the CP asymmetries are proportional to sin 2/3 (for B^ il^Kg, etc), sin 2a (for 
Bd —>■ TT+TT" etc.) or sin 2j (for Bg (pKg, etc). 

In the approximation that C2 — cs = 1 and that S1S2 ^ 1, we calculate sin 
2a, sin 2/3 and sin 27 to an accuracy of (1 +0(A^)), ie to 1% accuracy: 

sin2Q; = —2c^s^ 

sm2/3 = — -) 

si si 

sin27 =2c^s^4(l + ^)- (22) 

si C^Si 

It is interesting to note that S3 does not appear anywhere in these results. This 
is because all lengths of the unitarity triangle are simply proportional to S3. 
This is similar to the well known result that in the Wolfenstein approximation 
all lengths of the triangle are proportional to A. For us this lack of sensitivity 
to S3 is an essentially exact result. Given that we know Si precisely, and that 
is extracted from Si,S2 and the Cabibbo angle s^, the only uncertainties in 
numerically evaluating a, j3 and 7 come from experimental uncertainties in sin 9^ 
and in the dependence of S2 on 171^/171^, rric and 1]^ via x shown in equations 3 and 
4. We calculate sin 2a, sin 2/3 and sin 27 in terms of x for sin^,, = .221±.003. The 
results arc shown in the Figure. The solid line is for sin^^c = -221 while the long 
(short) dashed lines are for sin^c = .224(.218). Present experiments allow very 
wide ranges of a,/?, 7 : —1 < sin2Q;,sin27 < 1 and .1 < sin 2/3 < 1[19] so that 
our predictions are in a sufficiently narrow range that measurements of these 
CP asymmetries will provide a precision test of our model. Our predictions are 
very positive for experimentalists: sin 2/3 is not near its lower bound, and for the 
two most experimentally challenging cases, sin2Q; and sin 27, the asymmetries 
are close to being maximal. 
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How well can our model be tested with an asymmetric B factory operat- 
ing at the T(4S') with luminosity of 3.10^^ cm~'^s''^ [20]? We assume a total 
integrated luminosity of 10^^ cm~^, and find, using the numbers in [19], that 
for decay to a final state of branching ratio B, the quantity sin 2a (or sin 2/3 or 
sin 27) will be measured with an error bar ±S: 



where B = (4, 3,2)10"^ for Ba ij^Kg, B^ 7r+7r~,Ss pKg relevant for 
measuring sin 2/3, sin 2a and sin 27 respectively. Measuring all three quantities 
to ±.05 will provide a spectacular precision test of our model. The values of 
sin 2a, sin 2(5 and sin 27 must be fit by a single value of uiu/md which will be 
determined at the ±0.1 level. 

We stress two important features of our predictions for these CP asymmetry 
parameters. Firstly, as in the standard model, they are relatively insensitive to 
unknown QCD matrix elements. Secondly, they test the Georgi-Jarlskog ansatz 
in a deep way. For example the only dependence on the renormalization of 
gauge couplings beneath the grand scale comes from uncertainties in rjc- These 
uncertainties could be removed completely by taking the strange quark mass as 
input. Taking = 180 ± 60 MeV only leads to a 15% uncertainty in x- 

In this paper we have made accurate predictions for parameters in the 
neutral K and B systems, in the behef that the scheme of reference (1) will be 
decisively tested in the future. It is worth stressing that the predictions for B^ 
and VBJb are close to central quoted theoretical values, and thus are already 
strikingly successful. We have followed the consequences of the only framework 
incorporating the Georgi-Jarlskog mechanism which uses the minimal number 
of Yukawa couplings and has independent up and down quark mass matrices: 
there is absolutely no guarantee that e or Xd will be predicted correctly. Consider 
for example the case when e is dominated by the top quark contribution which 
is proportional to m^JBKRe{VtdV^lVusV*^). Even though a theory may give 
successful predictions for mt(100 - 200GeV), ^^^(.OOS - .018) and ^^.(.OSO - 
.054) it is not guaranteed that the prediction for e will be anywhere close to 
experiment. The quantity Re{VtdV^*VusV*^) has a spread of a factor of 50, 
and the quantity J which is proportional to SiS2s|s<^ could vary over a very wide 
range. In particular recall that the phase is determined by the requirement that 
\Vus\ = -221 ± .003. We think that it is extremely non-trivial that the prediction 
of our theory for m^Re{VtdV^lVusV*^)J is such that the central prediction for B^ 
is 0.4. 

The essential results of this paper are given in equation (16) for Bk (from 
e), equation (20) for VB/b (from Xd), equation (21) for Xg and in the Figure 




(23) 
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for sin 2ci;, 2/3, 27. The prediction for e'/e in equation (17) is less important as 
it involves uncertainties from the matrix elements. Once the top quark mass is 
accurately known, the range of predicted values for and VB/b will narrow. 
The largest uncertainly in comes from mu/md- CP asymmetries in decays 
of neutral B mesons offer the hope of a precision test of our theory which is free 
of strong interaction uncertainties. An asymmetric B factory operating at the 
T (4(5') with an integrated luminosity of 10^^ cm~^ can determine sin 2q;, 7 to 
an accuracy of ±0.05, and this will lead to a determination of mu/md to within 
±0.1. 
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